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We perform time-dependent simulations of spin exchange for an electron pair in laterally coupled 
quantum dots. The calculation is based on configuration interaction scheme accounting for spin- 
orbit (SO) coupling and electron-electron interaction in a numerically exact way. Noninteracting 
electrons exchange orientations of their spins in a manner that can be understood by interdot 
tunneling associated with spin precession in an effective SO magnetic field that results in anisotropy 
of the spin swap. The Coulomb interaction blocks the electron transfer between the dots but the spin 
transfer and spin precession due to SO coupling is still observed. The electron-electron interaction 
additionally induces an appearance of spin components in the direction of the effective SO magnetic 
field which are opposite in both dots. Simulations indicate that the isotropy of the spin swap is 
restored for equal Dresselhaus and Rashba constants and properly oriented dots. 
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I. INTRODUCTION 

One of ideas for solid state quantum computation em- 
ploys spins of electrons confined in quantum dotsi^ for 
information storage and processing. This idea drives an 
extensive experimental research on spin manipulation, 3 
detection^, decoherence and relaxation^ Construction 
of a universal quantum gate requires implementation of a 
controllable spin exchange operations between electrons 
confined in adjacent dots. In the absence of spin-orbit 
(SO) coupling, the spin-spin Hamiltonian possesses an 
isotropic form^ i.e. depends only on relative orientation 
of the spins. Anisotropic corrections^— are introduced 
by SO coupling. The spin processing procedures were 
revisec&i£ in order to remove or minimize the anisotropy 
due to the SO coupling. On the other hand a practical 
use was invented for the anisotropic exchange interaction 
in construction of quantum gates A recent studyi^ 
reported that the anisotropic part of the exchange inter- 
action vanishes in zero magnetic field which should lift 
the limitations for spin information processing that were 
the major concern of the previous workA^r— The con- 
clusion was supported by comparison^ of the double dot 
energy spectrum as found by the exact diagonalization 
technique and by a model Hamiltonian. 

The spin interactions for two-electron systems are 
probed by charging experiments^ that resolve the 
singlet-triplet avoided crossings due to the SO coupling. 
These avoided crossings occur in external magnetic field, 
for which anisotropy of the exchange interaction is evi- 
dent. At the moment the spin dynamics for double quan- 
tum dots in zero magnetic field can only be verified in a 
numerical experiment which we provide in this work. We 
present results of simulation in which the spin dynamics 
is monitored in time. We use the method of configuration 
interaction to simulate the spin swap in laterally coupled 
quantum dots. Our numerically exact results indicate 
that the swap process and result depend on the initial ori- 



entation of the spins also in zero magnetic field. We point 
out that the main source of the anisotropy is the effective 
magnetic field due to the SO couplin g 16 ! 17 that leads to 
precession of spins of moving electrons. The Coulomb 
interaction blocks the single-electron motion within the 
double dot. Nevertheless, the collective motion is still 
observed and we find that it results in the transfer of 
the spin associated with its precession. Moreover, we 
demonstrate that the Coulomb interaction leads to an 
appearance of spin components in the direction of the 
effective magnetic field which are opposite in both dots. 

The effects of the SO coupling for electron energy spec- 
tra is lifted and SU(2) symmetry is restored when the 
Dresselhaus and Rashba interactions acquire the same 
strength^ The Rashba interaction can be controlled by 
external electric fields and restoration of SU(2) symme- 
try allows for appearance of helical spin density waves 
in the two-dimensional electron gas^ Our simulations 
indicate that for equal Dresselhaus and Rashba coupling 
constants the spin swap becomes isotropic for carefully 
chosen spatial orientation of the double dot for which the 
SO effective magnetic field vanishes. 



II. THEORY 

We consider a two-dimensional Hamiltonian: 

h = + V(r)j 1 + Hsia + H BIA , (1) 

where k = — zV, 1 is the identity matrix, V(r) stands 
for the confinement potential, Hsia and Hbia intro- 
duce the Rashba and Dresselhaus SO interactions, re- 
spectively. For x and y axes oriented parallel to [100] 
and [010] crystal directions, SO interaction terms are 
given by H B ia = P (cr x k x - a y k y ) + H cub , and H S ia = 
^(c T x ky — a y k x )-\-H ( iiag- The cubic Dresselhaus term R cu h 
and the diagonal Rashba terms Hdiag (for their form see 
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|2Q| ) are included in the calculation but have a negligi- 
ble influence on the spin evolution. We use Ino.5Gao.5As 
effective mass m* = 0.0465mo, dielectric constant 
e = 13.55 and assume a model confinement potential 2 ^ 
V(x,y) = -V /[(l + [xVRl] 11 ) (1 + [y 2 /#*]")] + 

v b /[(i + [xVRl] 11 ) (1 + [yV^n where » = 10 > v ° = 

50 meV is depth of dots, 14 = 10 meV is the height of 
the inter dot barrier. The size of the double dot system 
in x and y direction is 2R X = 90 nm and 2R y = 40 nm, 
respectively. 2R b = 10 nm is the interdot barrier width. 

Calculations for the two-electron system start by diag- 
onalization of the Hamiltonian 



For each time step we evaluate the righthand side of Eq. 
(4) using the average values as provided by the time evo- 
lution in the basis of single-electron eigenstates 
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in the basis of eigenstates of Eq. (1) that are in turn de- 
termined using a multicenter Gaussian variational wave 
function^ The eigenvalues E m of Hamiltonian (J2j) and 
the corresponding eigenfunctions \I/ m (ri, cri, r2, a 2) ob- 
tained from the configuration interaction scheme are used 
for simulation of the time evolution as described by the 
Schrodinger equation ih^- = H^, 
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where c m coefficients are determined by the initial con- 
dition. For diagonalization of Hamiltonian (j2j) we use a 
basis of 325 two-electron wave functions obtained from 
26 lowest-energy single-electron eigenstates. The conver- 
gence of the results is discussed in the Appendix. 

In order to simulate the spin swap, in the ini- 
tial condition we localize the electrons in separate 
dots with opposite spin orientations. We denote the 
initial spatial single-electron wave functions localized 
in the left and right dots by ^ and ?/v, respec- 
tively. Functions and ip r are obtained by super- 
position of bonding and antibonding orbitals that are 
obtained as two lowest-energy eigenstates of Hamilto- 
nian (1) without the SO coupling. The initial two- 
electron wave function is taken as an antisymmetrized 
product = ^(^(ri)s + ((7i)^ r (r 2 )s_(a-2) - 

'0z(r2)5 + (<j2)V ; r(i , i)5_(cri)), where s + and s_ are orthog- 
onal eigenfunctions of a chosen spin component. Projec- 
tion of these wave functions on eigenfunctions of Hamil- 
tonian (2) defines c m = (^z^rtl^m) used in Eq. (3). 

Below we consider precession of a single electron spin 
in the effective magnetic field due to the SO coupling. 16 
For this purpose we solve the Bloch equation 



d(s) . . _ 

— =g/i b (B S o x s)/ft, 
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where is the Bohr magneton and ~Bso is the effective 
magnetic field due to the linear (dominating) terms of 
the SO coupling 
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FIG. 1: Spin swap simulation for quantum dots placed on 
the x axis. In the initial condition one electron is localized 
in the left dot (x < 0) with spin oriented antiparallel to the 
x axis and the other electron is localized in the right dot 
(x > 0) with the spin parallel to the x axis. Results were 
obtained for pure Dresselhaus coupling, (a) Black lines show 
the x component of the spin stored in the left (dashed) and 
right (solid) quantum dots.— The circles indicate the results 
obtained without SO coupling. Average values of y and z 
components is zero with (red lines) or without (red circles) 
SO coupling, (b) Spin-left (left column) and spin-right (right 
column) densities in selected moments in time. 



III. RESULTS 

Two-dimensional SO coupling constants depend on the 
thickness d of the confinement layer f3 = 7(7r/<i) 2 , and 
on the value of vertical electric field F z (external or built- 
in) a = qlsdFz, where 7 and asD are bulk Dresselhaus 
and Rashba constants, respectively. We use 7 = 32.2 
meVnm 3 for the InGaAs alloy 21 and assume d = 5.4 nm 
which gives f3 = 10.8 meV nm. This seems a maximal 
value of the coupling constant that can be practically 
achieved in a InGaAs quantum dot. The bulk Rashba 
constant of the alloy is aso = 0.572 nm 2 (after Ref. [22^. 
The electric field needed to produce a = 10.8 meV nm 
is 18.9 meV/nm, which is equivalent to the confinement 
potential drop of 102 meV across the dot of the height 
d = 5.4 nm. Results presented below stay qualitatively 
unaffected for smaller values of the coupling constants 
or weaker interdot coupling (for the latter - see the Ap- 
pendix). 
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FIG. 2: Simulation similar to the one presented in Fig. [T]only 
for electron spins set initially parallel or antiparallel to the z 
axis. In the left (right) column of plots the Coulomb interac- 
tion is included (neglected). (a,b) The average position of the 
spin- up density (black solid line). The electron-electron dis- 
tance (dashed lines) in the x (black), and y (blue) directions, 
calculated as \J {(x± — X2) 2 ) and ^ {(yi — 1/2) 2 ), respectively. 
(c,d) The spin components stored in the left (dashed curves) 
and right (solid curves) quantum dots.— Here and in all the 
other plots the x, y and z components of the spin are given 
by black, blue, and red lines, respectively. Circles in (d) - see 
text. 
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FIG. 3: Simulation for a single electron, pure Dresselhaus 
coupling and dots placed on the x axis. The electron is ini- 
tially localized in the left dot with spin oriented antiparallel 
to the x, y and z axis in the left, central and right columns of 
plots, respectively. (a,b,c) Average values of spin components 
are given by curves. The circles indicate the results obtained 
from the Bloch equation (4). (d,e,f) The spin stored in the 
left (dashed lines) and the right (solid lines) quantum dots. 
Black, blue and red lines show the results for x, y and z com- 
ponents, respectively. (g,h,i) The average x position of the 
electron packet. 



A. Quantum dots placed along [100] direction 

Let us first assume that the confinement potential is 
symmetric in the growth direction (a = 0) and that the 
centers of the dots are placed on the x axis. Initially the 
spin in the left (right) dot is set parallel (antiparallel) to 
the x axis. Figure Ufa) shows the time dependence of 
the average spin stored in the left and right dot^ For 
t = the spins in the left and right dots are not exactly 
equal to ±ft/2 due to leakage of ipi (ip r ) functions to the 
right (left) dot [see Fig. QJb) for t = 0]. The spin swap is 
completed at t = t s = 10.9 ps and no other component 
of the spin is generated during the process. 

The swap process exactly as illustrated in Fig. [1(a) is 
obtained in the absence of SO interaction independent of 
the choice of the direction in which the electron spins are 
initially set antiparallel to each other. In order to present 
the effects of SO coupling for the spin swap, let us change 
the initial condition. At t = the spin in the right (left) 
dot is now set parallel (antiparallel) to the z axis [see 
Fig. [2](c)]. At t = t s the absolute values of (s z ) are 
visibly reduced as compared to the initial condition and 
the spins in both dots acquire an identical non-zero value 
of the y component. Moreover, opposite components of 
the spin in the x direction are generated in both dots 
during the swap. The x spin components are maximal at 
t = t s /2 and disappear once the swap is completed. 

For a = the SO effective magnetic field of Eq. (5) 
is oriented along the x axis - the direction of electron 
tunneling between the dots. In order to evaluate the 
effects of the spin precession we performed calculations 
for a single electron. We localize the electron in the left 
dot and assume that its spin is oriented antiparallel to 



the x [Fig. [3(a,d,g)], y [Fig. [3(b,e,h)], and z axes [Fig. 
[3(c,f,i)]. The electron wave packet oscillates between the 
left and right dots [Fig. [3(g,h,i)]. The electron - localized 
initially in the left dot with the spin aligned with the z 
(y) direction - acquires a non-zero y (z) component of the 
spin when it tunnels to the right dot [Fig. [3(e) and [3(f)]. 
Figures [3( a, b,c) show that the spin evolution as obtained 
by integration of the Bloch equation (4) describing the 
spin precession perfectly agrees with the results of the 
main simulation. For the spin initially set antiparallel 
to the x axis - i.e. aligned with B50 - no precession is 
observed [Fig. [3(a,d)]. 

Let us now return to the problem of two electrons with 
spins initially parallel and antiparallel to the z-axis. The 
results for noninteracting pair of electrons [given by lines 
in Fig. [2(d)] are exactly reproduced by the sum of single- 
electron solutions [dots in Fig. Eld)]. Figure [2(b) shows 
that the horizontal electron-electron distance oscillates 
as the noninteracting electrons tunnel from one dot to 
the other and meet at the inter dot barrier. The same 
y component of the spin is generated in both dots in 
any moment of time. The electron spins are initially 
oppositely oriented with respect to the z axis and move 
in the opposite directions along the x axis. Accordingly, 
the change of y spin component as given by Eq. (6) has 
the same sign for both electrons. 

When the Coulomb repulsion is switched on, the oscil- 
lation of the electron-electron distance is no longer ob- 
served [Fig. [2(a)] - the electrons occupy fixed positions 
in separate dots. The electron density - the sum of spin- 
up and spin-down electron densities - is nearly station- 
ary, but the components of the sum are not. In Fig. 
[2(a,b) with the black solid line we plotted the center-of- 
mass of the spin- up electron density packet. We notice 



that this center oscillates in a very similar way for both 
interacting and noninter acting electrons. Also the spin 
swap as obtained for interacting electrons is similar to 
the one found in the absence of the Coulomb repulsion 
[cf. Fig. [2(c) and (d)], only the swap time is increased 
by a factor of 10 as the Coulomb repulsion enhances the 
effective height of the inter dot barrier. The only quali- 
tative feature introduced by the Coulomb interaction is 
the noticeable oscillation of the x component of the spin. 
We found as a general rule for interacting electrons that 
during the spin precession opposite spin components in 
the direction of ~Bso appear in both dots. 

Figure [4] shows the results for spins initially antiparal- 
lel in the y direction, still for the pure Dresselhaus cou- 
pling. The appearance of the z component of the spin 
during the swap is due to spin precession and is observed 
for both interacting [Fig. Ufa)] and noninteracting [Fig. 
ffl[b)] electrons. For the Coulomb interaction present the 
opposite spin components in the x direction (aligned with 
B50) appear in the dots, similarly as presented in Fig. 
[2(c) for the spins initially aligned with the z axis. 
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FIG. 4: (a) Same as Fig. [2(c). (b) Same as Fig. [2(d) only for 
electron spins initially antiparallel in the y direction. 

To further illustrate the findings of the above para- 
graph let us consider the case of pure Rashba coupling 
(Fig. [5]). In III-V material the Dresselhaus coupling can 
not be completely removed. Nevertheless, it can be small 
as compared to the Rashba coupling provided that the 
dots height is large and / or strong electric field is applied 
in the growth direction. For the pure Rashba coupling 
and the considered orientation of the dots ~Bso is aligned 
with the y axis [see Eq.(5)]. For the spins initially par- 
allel and antiparallel to this axis, the spin swap [Fig. 
E(c,d)] occurs without generation of neither x nor z spin 
components. For the spins initially aligned with x (z) 
axis one observes appearance of z (x) spin component - 
the same in both dots - for both interacting and non- 
interacting electrons - see Fig. [5(a,b) [Fig. E(e,f)] that 
results from the spin precession. According to Eq. (6) 
for f3 = the electron spin that is initially parallel to the 
x axis and moves in the x direction acquires the z com- 
ponent of the same sign as the spin antiparallel to the x 
axis that moves in the — x direction [Fig. [5(a,b)]. When 
the electron-electron interaction is present we addition- 
ally observe [see Figs. [5(a) and [5(e)] an appearance of 
opposite spin components in the y direction (i.e. direc- 
tion of B50) in both dots. 
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FIG. 5: Simulation for pure Rashba coupling and centers of 
dots placed on the x axis. The electrons in the left and right 
dots initially possess opposite spins in x (a,b), y (c,d), and 
z (e,f) directions. Black, blue and red curves show the x, 
y and z spin components stored in the left (dashed curve) 
and the right (solid curve) dots. Plots (b,d,f) were obtained 
for neglected electron-electron interaction which is included 
in (a,c,e). 
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FIG. 6: Simulation for both SO coupling types present a — 
ft/ 2 — 5.4 meV nm and centers of dots placed on the x di- 
rection. The electrons in the left and right dots initially pos- 
sess opposite spins along z ([001]) direction in (a,b), x/2 + y 
([1,2,0]) in (c,d), andx-y/2 ([2,1,0]) m (e,f). Purple, green 
and red curves show the [1,2,0], [2,1,0] and [0,0,1] spin 
components stored in the left (dashed curve) and the right 
(solid curve) dots. Plots (b,d,f) were obtained for neglected 
electron-electron interaction which is included in (a,c,e). 



Finally, Fig. [6] shows the results for both coupling 
types present with f3 = 10.8 meV nm and a = 5.4 meV 
nm. For the electrons moving between the two dots, the 
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SO magnetic field vector is aligned with the x — y/2 line, 
i.e. the [210] crystal direction. For the discussion of the 
spin behavior we consider this crystal direction and two 
orthogonal ones: within the plane of confinement [120], 
and in the growth direction [001]. When the electrons 
spins are initially aligned with the SO field vector, the 
swap goes without generation of other components of the 
spin [Fig. [6fe,f)] . For spins initially antiparallel to the 
[001] direction [Fig. [6(a,b)] the precession of the spins 
leads to an appearance of [120] spin component and vice 
versa [Fig. [6(c,d)]. The electron-electron interaction for 
both initial spin orientations that are orthogonal to B^o 
leads to appearance of the spin components parallel to 
the B^o direction [see the green curves in Fig. [6(a) and 
(e)] in the form noticed above for the pure Dresselhaus 
and Rashba coupling: vanishing at t = t s , maximal at 
t = t s /2, and always opposite in both dots. 

B. Spin oscillations due to the electron-electron 
interaction 

The results presented above indicate that during the 
swap of spins that are initially perpendicular to B^o, the 
electron-electron interaction induces appearance of spin 
components in the direction of this vector that are oppo- 
site in both dots [Fig. 2(c), Fig. 4(a), Fig. 5(a,e), Fig. 
6(a,b)]. The net spin in the direction of B^o remains 
zero, in contrast to the spin generated in the direction 
perpendicular to B50 by the spin precession. The spin 
components in ~Bso direction vanish at the end of the 
swap t = t 3 , however they are maximal at t = t s /2. 
Note that the XOR gate employs the square-root-of-a- 
swap operation^ i.e. the spin swap process interrupted 
exactly at t s /2. 

Let us analyze the background of the appearance of 
spin components in the direction of ~Bso vector. We 
focus on the simulation presented in Fig. 2(c), for which 
a = 0, the dots are placed along the x axis hence B50 
is aligned with the x axis, and the spins are initially set 
parallel and antiparallel to the z axis. 

By the Ehrenfest theorem the average electron-electron 
distance in the x direction changes in time as 

^-((X 1 -X 2 ) 2 ) = -^-(X 1 p xl -X 2 p x l) 

at m* 

4/3 2ih 

+ -T-{(T x lXl - (T x lX 2 ) -, (7) 

h m* 

where the last term compensates for the imaginary part 
due to the non-Hermitian xip x i operator. The electron- 
electron distance as presented in Fig. 2(a) is nearly 
constant but contains a rapid oscillation of small ampli- 
tude which results from a difference between the electron- 
electron separation in the initial condition and the equi- 
librium distance for interacting electrons (see Appendix) . 
For the purpose of analysis of Eq. (7) it is useful to limit 
the basis used for Eq. (3) to 4 lowest-energy two-electron 
states, which correctly describes the spin evolution (see 



Appendix) but is free of this rapid oscillation. The results 
for the spin evolution, electron-electron distance and the 
right-hand-side terms of Eq. (7) are displayed in Fig. 
For noninteracting electrons, at the right-hand side 
of Eq. (7) only the (x2P x i) = (x2)(p x i) term has a non- 
zero real part [Fig. [7(f)], which oscillates due to indepen- 
dent tunneling of electrons which go from one dot to the 
other with periodically changing positions and momenta. 
On the other hand the average value of x±p x i operator 
is purely imaginary. 24 For noninteracting electrons the 
term of Eq. (7) containing x<j x vanishes [Fig. 7(f)] and 
so does a x in both dots [Fig. 7(b)]. Thus the oscillation 
of the electron-electron distance observed in Fig. [7(d) is 
only due to the mean value of X2p x i operator. 

Interacting electrons keep their relative distance [Fig. 
[7(c)], so the terms at the righthand side of Eq. (7) must 
cancel one another. They indeed do [Fig. [7(e)] . Remark- 
ably, in contrast to the case without Coulomb interaction 
[Fig. [7(f)], for interacting electrons one finds [Fig. [7(e)] 

(xiPxi) = -(x2Pxi) (8) 

and 

(cr x ixi) = -(a xl x 2 ). (9) 

Relations (8) and (9) can be explained by analysis of 
the electron motion which becomes collective when the 
electron-electron interaction is introduced. The two- 
electron wave function of Eq. (3), can be written as a 
four-component wave function-^ V P 1 = ^^(ri, r 2 ), ^2 = 

tf n (n,r 2 ), ^3 = ^ t ( r i> r 2), and ^4 = ^u( r i> r 2), 
each corresponding to a given direction of the spin for 
a given electron label [e.g. \Ij^(ri,r 2 ) corresponds to 
electron of position ri (r 2 ) with spin oriented parallel 
(antiparallel) to the z axis]. Antisymmetry of the wave 
function with respect to the electron interchange implies 
*i(ri,r 2 ) = -^i(r 2 ,ri), tf 4 (ri,r 2 ) = -tf 4 (r 2 ,ri), and 
^2(1*1,1*2) = — ^3(r 2 ,ri). For the exchange of initially 
opposite spins ^3 and ^4 components are most relevant. 
Snapshots of \^s\ 2 and |^4| 2 are displayed in Fig. [8] as 
functions of x\ and x 2 calculated along the axis of the 
double dot 2/1 = 2/2 = for noninteracting [Fig. [8(a)] 
and interacting [Fig. [8(b)] electrons. Figure [8] shows 
that in the initial condition the electrons occupy sepa- 
rate dots and that the spin contained in the right (left) 
dot is oriented parallel (antiparallel) to the z direction. 
Spin orientation is inverted after the swap (t = t s ). For 
noninteracting electrons at t = t s /4 and t = t s /2 we no- 
tice [Fig. [8(a)] that probabilities to find both electrons in 
the same dot (i.e. on the diagonal x\ = x 2 of the plots) 
is non-zero, which results from an independent electron 
tunneling between the dots. Without the Coulomb in- 
teraction the spin swap occurs as due to single-electron 
tunneling. 

The interacting electrons do not occupy the same dot 
[see the vanishing probability density on the diagonal of 
plots presented in Fig. 8(b)] and the single-electron in- 
terdot tunneling is blocked. The interdot tunneling of 
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separate wave function components is still observed [Fig. 
8(b) for t = t s /4 an t = t s /2] but it occurs only along 
the antidiagonal of the plot X2 = —x\. Therefore, one 
can replace X2 by —x\ in the righthand side of Eqs. (8) 
and (9) which explains these relations. For the electron- 
electron distance to be unchanged the terms (xip x i) and 
(—X2Pxi) of Eq. (7) need to be canceled by the terms 
that contain the x component of the spin. The operator 
Xi(j Xl produces a non-zero contribution since the x spin 
component generated in the left dot (x < 0) has opposite 
sign [see Fig. 7(a)] than the one generated in the right 
dot (x > 0). We conclude that the generation of opposite 
spin components in the direction of the effective magnetic 
field is a consequence of fixed electron-electron distance 
and collective evolution of the two-electron wave function 
that are both induced by the Coulomb interaction. 
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FIG. 7: (a,b) same as Fig. [2c, d) but for the basis limited to 
four lowest-energy two-electron states. (c,d) Electron-electron 
distance (((xi — X2) 2 } 1 ^ 2 ) as calculated for the wave function 
given by Eq. (4) (black line) and integrated using Eq. (7) 
(red crosses) starting from the initial condition. (e,f) Real 
parts of righthand side terms of Eq. (7). 



C. Spin exchange for a — f3 

For a = f3 the linear SO terms of Hamiltonian (1) 
commute^ with -^(c x — &y) operator and the effective 

magnetic field B50 = §^(^x + k y )(l, — 1,0) T is aligned 
with [110] crystal direction. We performed simulations 
of the spin swap for a = ft = 10.8 meV nm. For the dots 
aligned in the x direction that were discussed above the 
results are qualitatively identical to the ones presented 
in Fig. 6 only with redefined direction of ~Bso vector. 
For a = () the direction of ~Bso does not depend on the 
orientation of the dots, however the strength of this field 
does. For the diagonal ([110]) orientation of the dots (the 
lower row of plots in Fig. [9J the electrons tunnel between 
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FIG. 8: Probability densities for the components of the two- 
electron wave functions corresponding to opposite spin orien- 
tations plotted on xi,X2 plane along the axis of the system 
2/1=2/2 = for chosen moments in time during the spin swap. 
Parameters are the same as in Fig. [71 The spins are initially 
oriented parallel and antiparallel to the z axis for noninter- 
acting (a) and interacting electrons (b). The spin swap time 
is t s = 1-48 ps for noninteracting electrons (a) and t s — 10.9 
ps for the interacting ones (b). 



the dots with equal k x and k y wave vectors, so the field 
should be relatively the strongest. On the other hand 
for the quantum dots oriented along the [110] direction 
(the upper row of plots in Fig. [9]) k x and k y have opposite 
sign for the electrons tunneling between the dots, so B50 
should be expected to vanish. 

The simulations show that spin exchange occurs in the 
same manner for the diagonal ([110]) and antidiagonal 
([110]) orientation of the dots only when the spins of 
electrons are initially antiparallel in the direction of B50 
([110] or x — y) - see Fig. 9(d) and 9(h). For the diagonal 
orientation of the dots and initial alignment of the spins 
in the direction perpendicular to the B50 vector: in z 
direction [Fig. 9(f)] and in [x + y) direction [Fig. 9(g)] 
one observes generation of (x+y) and z spin components, 
respectively. The precession is accompanied by genera- 
tion of opposite spin components in ~Bso direction [Figs. 
9(f,g)] in consistence with the results discussed above. 

Note, that the anisotropy of the swap for the diagonal 
orientation of the dots is observed in spite of the fact that 
for a = f3 the energy spectrum of SO-coupled system 
is identical 18 to the one obtained in the absence of SO 
coupling. For a = (3 the SO coupling does not affect the 
energy spectrum at zero magnetic field but the effective 
SO magnetic field is still present. 

For the antidiagonal ([110]) orientation of the dots (the 
upper row of plots in Fig. [9]) the spin swap becomes 
perfectly isotropic and occurs in the same manner inde- 
pendent of the initial spin orientation. Not a trace of 
spin precession is present in accordance with the single- 
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FIG. 9: Spin stored in separate dots for two interacting electrons with equal linear Dresselhaus and Rashba constants a — f5 — 
10.8 meV nm. The first and second row of plots correspond to different alignments of the double dot as shown schematically 
at the right end of the figure. In (a,e) (b,d) the spins are initially set parallel or antiparallel to the x, and z axes, respectively. 
In (c,g) the spins are initially set in the x + y (i.e. [110]) direction - within the plane of confinement and perpendicular to 
Bso direction ([110]). In (d,h) the spins are initially aligned with Bso vector. The dashed (solid) lines show the results for 
the left (right) dot. Black, blue, red, purple and green colors correspond to x, y, z, x + y and x — y components of the spin, 
respectively. 



electron picture of the electron tunneling that goes with 
k x — —k y in any moment in time which implies Bso — 0. 



D. Discussion 

The present study indicates that the spin swap as orig- 
inally defined for the purpose of controllable coupling of 
spin qubits^ localized in separate quantum dots is gener- 
ally anisotropic when the spin-orbit coupling is present. 
The anisotropy of the spin swap results from the effective 
magnetic field due to the spin-orbit coupling. This field 
changes the direction of the electron spin as it moves 
in space. The study of Ref. [3 indicated that for a 
carefully chosen computational two-electron basis the ex- 
change Hamiltonian becomes formally isotropic at zero 
magnetic field. The proposed^ basis is obtained by a 
unitary transformation of a separable basis of singlet and 
triplet states. The unitary transformation [Eq. (13) of 
Ref. 14] produces basis elements in which the spin and 
spatial degrees of freedom are entangled, i.e. direction of 
the electron spin depends on its position in space. For 
the purpose of the quantum computation any basis can 
in principle be chosen. However, the entangled basis that 
allows for a simpler form of the Hamiltonian requires a 
more challenging handling of the quantum information, 
which in fact should be stored by entangled spin-orbital 
wave functions rather than by the electron spin itself. 
The practical usage of the entangled basis calls for new 
physical procedures for preparation of the initial state 
and read out of the quantum computation result. 

For construction of the universal quantum gate the 
two-spin operations need to be complemented by single 
spin rotations. The direct idea to perform the latter is 
to put the system in external magnetic field to split the 
spin states of the single-electron and exploit the Rabi 



oscillations in resonant radiation of microwave or radio 
frequency^ In presence of the external magnetic field 
(B) the spin swap becomes anisotropic even without SO 
coupling, since the electron spins precess in B unless they 
are initially aligned with the direction of the external 
field. In order not to interfere with the spin exchange 
the single-spin rotations should be applied without the 
external magnetic field. The original idea for that pur- 
pose was the electrically controlled coupling of a selected 
spin to a ferromagnetic medium^- It was also demon- 
strated that the single-spin rotations can be performed 
using the spin precession in the SO effective magnetic 
field which occurs when the electron is made to move, 
e.g. along closed trajectories^ This idea for the single- 
spin rotations^ does not require the external magnetic 
field or irradiation. 



IV. SUMMARY AND CONCLUSIONS 

We presented numerically exact simulations of the spin 
swap for two-electron SO-coupled double quantum dots. 
The study covered Dresselhaus and Rashba interactions 
and various spatial orientations of the double dot. The 
swap of spins as observed in time-dependent simulation 
involves four mechanisms: (i) direct tunneling which con- 
sists in electron carrying its spin from one dot to the 
other (ii) the spin tunneling which still occurs when the 
direct tunneling is blocked by the electron-electron inter- 
action (iii) the precession of the spin moving in the ef- 
fective magnetic field due to the spin-orbit coupling and 
(iv) generation of opposite spin components in the direc- 
tion of the effective magnetic field which is observed for 
interacting pair of electrons. The third and fourth mech- 
anisms of the above list can be switched off for initial 
spin orientation aligned with the effective magnetic field 
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vector B^o- For the initial orientation of the spins in 
one of the directions perpendicular to B50 the spin in 
the other perpendicular direction is generated during the 
swap as a consequence of the spin precession. We argued 
that mechanism (iv) is necessary to maintain a constant 
electron-electron distance and is a consequence of a col- 
lective motion of the electrons within the inner degrees of 
freedom which is still observed when the single-electron 
tunneling between the dots is blocked by the Coulomb 
repulsion. We also demonstrated that for both coupling 
types present the spin swap is largely affected by a spe- 
cific orientation of the double dot system within the (001) 
plane of confinement via the strength of the effective mag- 
netic field. In particular, we demonstrated that the SO 
coupling effects can be eliminated from the spin swap 
process for quantum dots aligned with [110] crystal di- 
rection. For this orientation of the dots and the Rashba 
constant tuned to match the Dresselhaus constant B50 
vanishes and the spin swap becomes perfectly isotropic. 
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FIG. 10: Results of simulations for the pure Dresselhaus cou- 
pling and dots placed on the x axis for the interacting elec- 
trons with spins initially oriented antiparallel in the z direc- 
tion. Black, blue and red curves show the spin components 
stored in the left (dashed curve) and the right (solid curve) 
dots. The green curve at the top of the plot is the average 
electron-electron distance that is referred to the right axis. 
The results presented in (a), (b), (c) and (d) correspond to the 
two-electron basis containing m = 4, 10, 50 and 325 lowest- 
energy eigenstates of Hamiltonian (2). 



Appendix 

This Appendix presents convergence of the results with 
respect to the energy cut-off used in Eq. (3). Let us con- 
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FIG. 11: Same as Fig. [TuT d) only for interdot barrier height 
increased from Vj, — 10 meV to Vj, — 50 meV. 



sider the dots placed on the x axis and electron spins 
initially oriented parallel and antiparallel to the z axis 
that was discussed in the context of Fig. 2(c). Fig. 
10 shows the results for m = 4, 10, 50 and 325 two- 
electron lowest-energy eigenstates used as basis elements 
in Eq. (3). The basis with m = 4 covers the ground-state 
and three-fold degenerate excited stated - which in the 
absence of SO coupling corresponds to the spin triplet. 
The energy separation of the ground-state and the ex- 
cited state is AE = 0.189 meV, which well corresponds 
to the spin swap time 2 t s = tt/AE = 10.9 ps. The ba- 
sis of 10, 50 and 325 elements covers all the two-electron 
eigenstates of the energy that exceeds the ground-state 
energy by not more than 8.9, 19.4 and 65.2 meV, respec- 
tively. For m > 4 a rapid and low amplitude oscillation 
appears in the results of Fig. 10. This oscillation results 
from a difference between the electron-electron separa- 
tion in the initial condition and the equilibrium distance 
for interacting electrons. The initial condition is taken 
from single-electron wave function obtained for nonin- 
teracting electrons which are projected onto the basis of 
two-electron eigenstates (see Section II). The electron- 
electron distance is nearly the same in the lowest-energy 
four states (singlet and triplet states). The constant 
electron-electron distance obtained for m = 4 [Fig. 10(a)] 
is the equilibrium distance for interacting electrons in the 
ground- state. For larger m the basis resolves between 
the equilibrium distances for interacting and noninter- 
acting electrons, hence the appearance of the rapid os- 
cillations of electron-electron distance and the resulting 
oscillations of the spins. The oscillations do not affect 
the mechanism of the spin swap nor the swap time and 
have a small amplitude which can be further reduced 
by a choice of confinement parameters. In particular 
Fig. 11 shows the results for the barrier height increased 
from 10 do 50 meV^ The Coulomb interaction affects 
weakly the electron-electron equilibrium distance and the 
rapid oscillations have a negligibly small amplitude also 
for m = 325. The results presented in this paper were 
obtained for the fully convergent 325 element basis with 
the exception of Subsection III.B, where we use the four- 
element basis for simplicity. For AE = 65.2 meV the 
shortest oscillation period that can be accounted for is 
0.06 ps. 
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